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DEVELOPMENTS IN COMPUTATIONAL MODELING OF 
TURBULENT FLOWS 
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lence will be presented. Interested readers may mathematically by using the generalized Cayley- 

reference the recent papers by Perot and Moin Hamilton formulas and the invariant principles 

[25, 26] which shed some physical insight on the [27]. In what follows we will, based on this general 

behavior of near-wall turbulence from their DNS constitutive relation, discuss the development of 

of shear-free turbulent boundary layer flows. Reynolds stress algebraic equation models. 
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where ever, in practice, one may not need all of these 

terms. For example, the quadratic tensorial form 
IIi = Ui, k U k ,i , II 2 = U i<k U i>k , j s quite appropriate for turbulent shear flows [13]. 

n 3 = Ui, k Ul k , n 4 = UlkUlk , Launder [31] found that the cubic terms in Eq. 

rb = Ui kU’ihUii , n 6 = UiuU? k Uh , (7) are needed for turbulent flows with swirling. 
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However, determination of all the coefficients is 2.2.1 or by the new k-e equations described in 
not an easy task. For practical purposes, Shih Section 2.2.2. 

et al [13] suggested a quadratic stress-strain re- The model represented in Eqs. (19)-(23) is 

lation, which is also suggested by an RNG result quite simple but has several advantages compared 

[29]. The necessity of using a higher order stress- to the standard k-e eddy viscosity model. First, 
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The O 2 and its dissipation rate eg will be deter- be written as 
lined by the following transport equations: -q 
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terms which are considered as the most impor- In addition, the modeling of the slow pressure- 

tant model terms in the second moment equa- strain correlation is traditionally combined with 

tions. Many attractive models have been pro- the dissipation rate tensor and the combination 
posed by using a variety of theories and principles, is called the return-to-isotropy term: 
such as invariant theory, rapid distortion theory 1 • 2 


• +9 4-9 

S 2 2 

rj 4-9 +9 

o <2 ” cl 

i s 

■sag 

a s-t* 


. ce 

X3 N 


C3 <L3 4^ 

0> fn CO 
TO 2 

2 8 S 

- | 8' 

”3 S 8 

49 

.2 >3 " 
■S 3 *3 
.gag 

/o. ju o 

f— 1 ci) 

. h H 4 ) 


03 

£ + 


m A 
TO S 

O 


T 

0 3 

K« 

*? 15 

H f-i 

! £ 


<3 Pi ^ 


LO 


N w 

C) 
f-t 

<D XJ 

■** •§ S 

111 

03 >i 

® -3 0 
81 § 
^ H ’•§ 

•a a | 

y S h 

3 ”® J§ 
£ ,0 +» 


TO ,H 


XJ 
rt p 

II 

.2 § 


■8 §• 


p S 

C7* 0 
0 u 




•©* ^ 
+ ^ 


cs ft. 

<N R, rO 
►C V ' 


, CW •<-, 
. *© w a 

54 -o 

R< •*» 
►O M V 
v ✓ rO 


co qj a> 
.SJ ,d 

8 ■* g 

s * 2 ' 

S.tf.3. 

$» ; 

g j* 

CO - 5 

Q> -2* ’3 
^ g ft 

X) ^ X> 


^ > 3 

►« p o 

*© I®* 43 

»-C> _ +i 

* hH 

II ® . 

X> 

c^ 4T ,<? <p 

*c> g, co 




ffl *s i? 

O 4i ^ 
O co O 


t_T io 

s as 
5 a . 

r> <-Q d< 


*o 4 =; 
H l<N cpj 


cu co 4-5 
Ph *43 
V 03 O 

'U W 45 


<0 

1 

kcT 

R. 

rO 


1 

rO 


>3 

cs-S* 

+ 

04 

04 

« ft. 
►O 

+ 

04 

i-H |CO 

4- 

04 

►0 

<<T 

1 

*cT 

+ 


+ ft 

csTTco 

N-S* + 


O X> 

^ 15 13 

•a £ 1 

p +> ST 


Ago 

,2*g a 

sSs 

8 h 42 

d is o> 

CO *3 f-4 

a I % 

rH -P 

•g § 1 

42 P< *t2 


S3 d p 

K <v m , 

l ? 3 a *8 

*8 50 5 i 
rt g| 

O O X5 

ts g>’8 

fl .3 o 

P > cx 


42 a or 

d r ^ 

1 H <2 


KgH 
« 6 


o S 

g & 

■4-5 R 

rW "g 

■ 9 I 

1 «-• d 


tH Q) 

co d > 

co ’43 *43 

xJ d 
^ 15 43 
o d 
d o> -g 


»H rr-J O 

«2 S « 

12 13 <v 

S > & 


I £ 
*3 

OJ o 

■$ a 


c6 d <d 

■3 .§ * 

m W c 3 
<1> 2 CO 
.6 

-4-5 O, 

« « £ 

ia«§ 

o 15 +=> 

o 3 d 

, rd <U 

<3 +5 4d 


■JL w f 

.2; a> g 

\ d o 

^li 
8 I ^ 

3 0 

3 u rl 

1 -a | 

s ! a 

" 4-5 ,-d 

! ^ *y 3 

, M rd O 

I £ w 


g O § 

M § 

HH ^ 45- 

co O co 

• t-l * ^ 


8 1 S 

.S *> g 

t> +5 4^ 


eo 

co 



3.3.1 Expansion in terms of I I b and III b much less than thirty two. Using the realizabil- 

As mentioned previously, the seven undeter- ity constraint Eq. (49), the number of the con- 

mined coefficients in Eq. (53) are, in general, stants can be substantially reduced. Johansson 

functions of Ilf, and IIIb- According to the defini- and Hallback found that only four undetermined 
tions of Eq. (54), the absolute values of II b and constants remain unresolved. The rapid pressure- 
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Fig. 14. Center velocity distribution in a diffuser 

r^2] : present model; : SKE; •: experiment. 









Fig. 17. Evolution of normalized temperature Fig. 19. Mean temperature profile for the con- 

variance and its dissipation rate (Sirivat and stant temperature flat plate boundary layer (Gib 

Warhaft [44]). son [45]). 
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Fig. 26. Direct comparison of the return-to-isotropy models with 
DNS data of homogeneous shear flow C128W (Rogers et al. [50]) 
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11. SUPPLEMENTARY NOTES 

Prepared for the International Symposium on Mathematical Modelling of Turbulent Flows sponsored by the University of 
Tokyo, Tokyo, Japan, December 18-20, 1995. ICOMP Program Director, Louis A. Povinelli, organization code 2600, 
(216) 433-5818. 
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